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Abstract. The temperature-induced second-order phase transition from 
Bloch to linear (Ising-like) domain walls in uniaxial ferromagnets is investi- 
gated for the model of D-component classical spin vectors in the limit D — » oo. 
This exactly soluble model is equivalent to the standard spherical model in the 
homogeneous case, but deviates from it and is free from unphysical behavior in 
a general inhomogeneous situation. It is shown that the thermal fluctuations 
of the transverse magnetization in the wall (the Bloch-wall order parameter) 
result in the diminishing of the wall transition temperature Tb in compari- 
son to its mean-field value, thus favouring the existence of linear walls. For 
finite values of Tb an additional anisotropy in the basis plane x,y is required; 
in purely uniaxial ferromagnets a domain wall behaves like a 2-dimensional 
system with a continuous spin symmetry and does not order into the Bloch 



1. Introduction 

The spherical model proposed in 1952 by Berlin and Kac (see also ||) is ex- 
tensively used up to now as the only exactly soluble model describing the phase 
transition in 3-dimensional magnetic systems. In contrast to the mean field approx- 
imation (MFA), the spherical model describes, in a simplified manner, the thermal 
fluctuations of spins, which can be taken into account exactly due to their Gaussian 
nature. The technique for the consideration of inhomogeneous systems described by 
the spherical model was developed by Barber and Fisher Q for the surface-induced 
inhomogeneity in layered magnetics and elaborated by Abraham and Robert || for 
the problem of phase separation (i.e., the domain wall (DW) formation) in the spher- 
ical model. Later the inhomogeneous states of the bounded spherical model induced 
by antiperiodic M and twisted |?J boundary conditions were investigated. 

The results obtained for inhomogeneous states of the spherical model posess some 
unexpected features. The phase transition temperature T c of a 4-dimensional ferro- 
magnetic slab consisting of N 3> 1 layers turns out to be higher than that of the bulk 
in the case of free-edge boundary conditions [Q. The two-domain state induced by 
the magnetic field ±H in two half-spaces is characterized by the DW width diverging 
and the domain magnetization vanishing in the limit H — > 0, i.e., in contrast to the 
underlying Ising model the spherical model does not exhibit the phase separation || . 
As was argued already in such features are the consequence of the global spin 
constraint |l| , which obviously becomes not so good in the inhomogeneous case. 
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Another version of the spherical model — the model of isotropic D-component 
classical spin vectors in the limit D — > oo — was proposed by Stanley [|, who 
showed that in the homogeneous case this model is equivalent to the spherical model 
by Berlin and Kac Q. The normalization condition for a spin on a lattice site i, |mj| = 
1, becomes, in the case D — > oo, very similar to the global spin constraint, which is the 
reason of the equivalence of the two models. However, since the spin normalization 
condition does not mix the spin variables on different lattice sites, the D- vector model 
by Stanley with D — > oo is more physically appealing than the original spherical 
model. Moreover, the two models become non-equivalent in a general inhomogeneous 
situation, where, as was shown by Knops j[o|, the D = oo model corresponds to some 
generalization of the spherical model using a local spin constraint. This idea was 
substantiated in the work by Costache et al fll|| , who calculated the Curie temperature 
T C (N) of a ferromagnetic film of N 3> 1 layers with free boundary conditions in 4 
dimensions using a set of independent spin constraints in each layer. The numerically 
calculated values of T C (N) monotonically increase with N to the bulk value TJ oo), 
which is physically expected and differs from the result of Barber and Fisher jf[ for 
the standard spherical model. 

The further advantages of the Z)-vector model are the possibilities of consideration 
of finite-D and anisotropic systems. The latter is important, in particular, for the 
calculation of finite-size corrections to T c of ferromagnetic films mentioned above in 
the actual 3-dimensional case. Since such a film with N < oo is a 2-dimensional 
system, the T c -corrections are finite only in the presence of the stabilizing uniaxial 
anisotropy ]l2| . In spite of its advantages in comparison to the standard spherical 
model, the D-vector model with D 3> 1 was much less used than mentioned. As 
exclusions one can cite the works by Abe and Hikami |l3|, [l4| and by Okabe and 
Masutani |l5| dealing with the 1/D expansion for 3-dimensional systems and the 
work by Okamoto Jl6| where the uniaxial spherical model with a transverse field was 
considered. 

It should be noted that practically all the researchers dealing with the isotropic 
D- vector model used the designations N or n instead of the original D. Such a 
modification is, however, not justified in a general anisotropic case, where the number 
n of the relevant order parameter components determining the symmetry and thus the 
universality class of a system can be different from D. As an example one can consider 
the rather general "n-D model" Jl7|, |f8) having the first n < D components coupled 
by the exchange interaction with the equal strength and the rest D — n components 
"free". Among realizations of the n-D model are, in particular, the x-y model (D = 3, 
n = 2) and the plane rotator one (n = D = 2) belonging to the same universality 
class determined by n but having different values of T c depending on both n and D. 
Correspondingly, in a general case the 1 jn expansion of the critical indices is not the 
same as the 1/D expansion of non-universal quantities. 

The general qualitative result of Q , the absence of the phase separation in the 
spherical model (but not the disappearance of the domain magnetization !), can be 
explaned by the fact that this model behaves in the bulk like the isotropic D = oo 
model |^|, [l?]], which obviously exhibits no phase separation. For the D- vector model 
the separation of a specimen into domains with opposite magnetizations by domain 
walls of a finite width requires an easy-axis anisotropy, which makes the intermediate 
orientation of the magnetization in the wall energetically unfavourable in comparison 
to that in domains. Clearly, this actual situation cannot be treated either with the 
help of the spherical model in its standard formulation J5J], or with the improved 
version |Tl], [T^], which is equivalent to the isotropic D — oo model in the general 
inhomogeneous case. 

The problem arising here — the study of the influence of thermal fluctuations 
on the domain wall structure — is not only important for comparing the properties 
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of different versions of the spherical model. The physics of domain walls at elevated 
temperatures is itself rather interesting and unexplored area, whereas since the sem- 
inal work by Landau and Lifshitz pOfl the majority of researchers have addressed 
the zero-temperature statics and dynamics of DWs based on the assumption of the 
constant magnitude of the magnetization in the wall. 

The first theoretical investigation of the temperature variation of the structure 
of domain walls is due to Bulaevskii and Ginzburg ]2l] j. who with the help of the 
phenomenological version of MFA, using a macroscopic Landau free energy in the 
vicinity of T c , predicted a phase transition from Bloch to linear DWs in uniaxial 
ferromagnets at some Tg < T c . Qualitatively this phase transition can be explaned 
by the fact that the spins in the center of a Bloch wall, which are forced to lie 
perpendicular to the easy axis, experience a molecular field smaller than in domains 
and hence order at some temperature Tg less than T c , which leads to linear (Ising-like) 
walls (LW) in the region Tg < T < T c . For ferromagnets whose anisotropy energy is 
much smaller than the exchange interaction, the LW temperature region is narrow. 

Complementary, the transition from Bloch to linear walls at T — depending on 
the anisotropy was investigated by van den Broek and Zijlstra |22]| with numerical 
methods. It was found that LWs are realized if the ratio of the anisotropy energy 
to the exchange one exceeds 2/3; the DW width S, is in this case, comparable with 
the lattice spacing ao- Later this transition was discovered by Sarker, Trullinger and 
Bishop p3| in the framework of a formal soliton theory independently of Bulaevskii 
and Ginzburg. The problem was addressed also by Niez and Lajzerowicz [[24|, p5[ , 
where the factor 2/3 mentioned above was calculated analytically. 

The first indirect experimental evidence for the transition from Bloch to linear do- 
main walls was obtained from the optical observations of the temperature dependence 
of the period of the domain structure in YFe03 just below T c ^6|. Lately the LWs 
were observed in the dynamical susceptibility experiments on the low-temperature 
ferromagnets GdCl 3 @ and LiTbF 4 @. In @ §(J the DW mobility was calcu- 
lated in the whole temperature range, which exibited a deep minimum at Tg. Such 
a minimum was observed recently in the dynamical susceptibility experiments on the 
high-temperature Ba and Sr hexaferrites ju| p2| . 

Recent experiments also provided evidence of strong fluctuational effects about 
the DW phase transition. The transition temperature Tg was substantially lower than 
its mean-field estimate, and the critical index j3g of the Bloch-wall order parameter 
(the transverse magnetization in the centre of the wall) was about 0.1 in contrast 
to the MFA value 1/2. Such strong fluctuations are actually not surprising since 
a domain wall is a 2-dimensional object. The analysis by Lawrie and Lowe (3^] 
making use of renormalization-group arguments has led to the clear result that a 
domain wall in a biaxial ferromagnetic model having an additional anisotropy in the 
x,y plane (which is usually the effective one due to the magnetostatic field 20 1) 
belongs to the universality class of the 2-dimensional Ising model, and hence one can 
expect (3g — 1/8. In contrast, in a purely uniaxial model without the dipole-dipole 
interaction a domain wall behaves like a 2-dimensional plane rotator model and can 
show only the Kosterlitz-Thouless phase transition without ordering to a Bloch wall. 

The absense of the long-range order (i.e., the transverse magnetization compo- 
nent) in a domain wall in a purely uniaxial ferromagnet can be demonstrated with 
the help of the linear spin- wave theory. The thermal disordering of Bloch walls is due 
to the so-called Winter magnons [ p4| , the excitations localized on the domain wall 
with the dispersion law e 2 oc Aq 2 {Acp + K±) (A is the inhomogeneous exchange con- 
stant, K± is the in-plane anisotropy constant and q is a 2-dimensional wavevector). 
The first factor in e 2 corresponds to the free translational motion of the wall and the 
second describes the rotation of the magnetization in the center of the wall in the x, y 
plane. It can be seen that the number of Winter magnons, which in the classical case 
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is proportional to J d 2 q/e q , diverges logarithmically at small q, if K± tends to zero. 
Thus, in a purely uniaxial ferromagnetic model the linear walls cannot order to the 
Bloch ones at any non-zero temperature. For a small non-zero in-plane anisotropy 
K± the shift of Tb from its MFA value due to fluctuations should be very essential. 

The aim of this work is to find an exact solution for the domain wall magneti- 
zation profile in biaxial ferromagnets at non-zero temperatures and their transition 
temperature Tb in the framework of the spherical model in its D- vector version. In- 
stead of modifying the approach of |B| based on the calculation of the partition 
function with the steepest descent method, we will use the diagram technique for 
classical spin systems [17| |3?|. This diagram technique, which is a generalization of 
the "Ising part" of the spin operator diagram technique of Vaks et al [^6[ [37J , makes 
it possible to locate and sum up all the diagrams surviving in the limit D — > oo and 
can be reformulated for the present purposes for inhomogeneous situations. 

The approximation obtained by summing up such diagrams (without going to 
the limit D — > oo) is the so-called self-consistent Gaussian approximation (SCGA), 
which was first formulated by Horwitz and Callen |}8| for the Ising model (D = 1). 
SCGA yields rather good results for the thermodynamic quantities of the Ising |3S[| 
and the classical Heisenberg (D = 3) |Q models on 3-dimensional lattices in the 
whole temperature range and can be of importance for a possible improvement of the 
presently obtained results for the domain wall structure in the spherical limit with 
regard to systems with finite D. A detailed description of the classical spin diagram 
technique and SCGA can be also found in the recent publication Q. 

The rest of the paper is organized as follows. In Section 2 the diagram technique 
for classical spin systems and the construction of SCGA in the inhomogeneous case 
is described. In Section 3 SCGA is simplified for D — > oo, and a closed system of 
equations for magnetization and spin-spin correlation function describing the domain 
wall in the spherical limit is derived. In Section 4 the magnetization profile of a fluc- 
tuating domain wall is calculated and the dependence of the transition temperature 
Tb on the anisotropy parameters is analyzed. In Section 5 further problems of the 
DW statics and dynamics at elevated temperatures are discussed. 



2. Classical spin diagram technique and SCGA 

The appropriate classical D-vector Hamiltonian with biaxially anisotropic ferromag- 
netic exchange interaction can be written in the form 

H = -- ^2 J iJ ( m zi m zj + ir m yi m yj + X! r \am ai m a3 J , (2.1) 

ij \ a=3 / 

where i,j are the lattice sites, is the normalized D-component vector, |rrii| = 1, 
the dimensionless anisotropy factors satisfy rj a < rj < 1 and all rj a are, for simplicity, 
taken to be equal to each other. For D = 3 equation (24) reduces to the anisotropic 
classical Heisenberg model. In the chosen geometry the average magnetization in the 
bulk is directed parallel or antiparallel to the easy axis z, in the centre of a Bloch wall 
it takes on one of two possible orientations along the second easy axis y. All variables 
describing the magnetization profile of a plane DW are functions of the coordinate x 
only. The temperature-normalized molecular field ^ i acting on a spin on the site i 
from its neighbours is given by 

8(BH) r ( ^ \ 

' = ~ c)m = ^ J i \ mzjez + r f m vj e y + Vam aj e a I , (2.2) 

1 j \ a=3 / 

where (3 = l/T and e are unit vectors in appropriate directions. The mean field 
approximation consists in neglecting fluctuations of the molecular field replacing 
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m z j =>■ (m z j), rriyj =>■ (m y j) and m a j =4> in ( |2.2| ), one arrives at the inhomogeneous 
Curie- Weiss equation 

(m. t ) = (2.3) 

where i?(£) is the Langevin function. For small anisotropy (1 — r\ <C 1) the magnetiza- 
tion varies slowly on the scale of the lattice spacing, and the continuous approximation 
can be applied to ( |2.3| ) . In this case the zero-temperature results of Landau and Lif- 
shitz [^0) and the finite-temperature ones of Bulaevskii and Ginzburg [ElJ for the DW 
magnetization profile are recovered (see below). 

Fluctuations of the molecular field (|]^) can be taken into account within the 
framework of a perturbative scheme based on the diagram technique for classical 
spin systems |L7|, 18 . The perturbative expansion of the thermal average of any 



quantity A characterizing a classical spin system (e.g., A — m,) can be obtained by 



rewriting (2.1) as TL = 7io + Vi n t, where Tio is the MFA Hamiltonian with the averaged 



molecular field determined by (|2.3|), and expanding the expression 



JV 



M = \ I II dmjAeM-PH), |m, | = 1 (2.4) 

in powers of Vi n t- The averages of various spin vector components a, /3, 7, ... = 
1, . . . , D on various lattice sites i,j,k,... with the Hamiltonian TIq can be expressed 
through spin cumulants, A .., (see below) in the following way: 

(m a i)o = A Q j, 

(m a imi3j)o = A a pi5ij + A m Apj, (2.5) 
{m ai mpjm 1 k)o — A a p 7i 5ijk + A a p i A 1 ]~5ij 

+ ApjjA ai 8jk + A iai A/3j5ki + A ai ApjA 1 k, 

etc., where 5ij, Sijk, etc., are the site Kronecker symbols equal to 1 for all site indices 
coinciding with each other and to zero in all other cases. For one-site averages (i — 
j = fc = . . .) equation ( |2.5| ) reduces to the well-known representation of moments 
through cumulants (semi-invariants), generalized for the multi-component case. In 
the graphical language (see, e.g., Fig. 0) the decomposition (2J5) corresponds to 



all possible groupings of small circles (spin components) into oval blocks (cumulant 
averages). The circles coming from Vint (the "inner" circles) are connected pairwise by 
the wavy interaction lines representing the quantity rj a /3Jij. In diagram expressions 
the summation over the coordinates i and component indices a of inner circles is 
carried out. One should not take into account disconnected diagrams (i.e., those 



containing disconnected parts with no "outer" circles belonging to A in (2.4)), since 
these diagrams are compensated for by the expansion of the partition function 2 in 
the denominator of ( |2.4D . Consideration of numerical factors shows that each diagram 
contains the factor l/n s , where n s is the number of the symmetry group elements of a 
diagram (see ( |2.9| ) and ( |2.10 ), the symmetry operations do not concern outer circles). 



In the homogeneous case it is more convenient for practical calculations to use the 
Fourier representation and calculate integrals over the Brillouin zone rather than 
lattice sums. As the lattice sums are subject to the constraint that the coordinates of 
the circles b elon ging to the same block coincide with each other (due to the Kronecker 



symbols in (2^5)), in the Fourier representation the sum of wavevectors coming to or 



goin g out of any block along interaction lines is zero. The cumulant spin averages in 



(2.5) can be obtained by differentiating the generating function A(£) = lnZo(£) over 



appropriate components of the dimensionless molecular field £: 



A (£) ^ (2 6) 
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Fig. 1: Self-consistent Gaussian approximation (SCGA) for classical spin systems, (a) 
and (c): block summations for the renormalized magnetization and pair spin cumulant 
averages; (b): Dyson equation for the spin-spin correlation function. 



where -Zo(£) = const x £~( d ' 2 ~ 1 'Id/2-i(Q is the partition function of a D-component 
classical spin and is the modified Bessel function. The two lowest-order cumu- 
lants which will be needed below can be written explicitly as 

A a (£) = m)j, (2-7) 

A Q /3(C) = —J— I a /3 12~ j + B (0-^-: 

where 5 a p is the spin component Kronecker symbol, D(£) = c?A(£)/d£ is the Langevin 
function for D-component spins and B'(£) = dB(£)/d£. The expressions for the 
3- and 4-spin cumulants can be found in |p5[. It should be stressed that the spin 



cumulants (2.7) appearing in the diagrams generated originally by the expansion of 
(2.4) in powers of Vi n t (the unrenormalized diargams) simplify, since there is only 
few non-zero components of the molecular field £ (for a domain wall in the chosen 
geometry £ z and £ y ). The complete form of spin cumulants fl2.7|) is needed, however, 
for the construction of SCGA allowing for the fluctuations of other components of 
the molecular field. For Ising systems the classical spin diagram technique coincides 



with the "Ising part" of the standard SDT |5G, 37 and can be used with Brillouin 
functions B$ of a general spin S. In the book [37| more technical details concerning 
the construction of SDT for Ising systems can be found, which play the same role in 
the present classical SDT. 

The next step beyond MFA is the self-consistent Gaussian approximation taking 
into account pair correlations of the molecular field acting on a given spin from its 
neighbors, which implies the Gaussian statistics of the molecular field fluctuations (see 
Fig. |l|). Since subsequently we are going to take the limit D — > oo, only fluctuations 
of the molecular field components with a = 3, . . . , D should be taken into account, 
because their total contribution exceeds that of the fluctuations of z and y components 
by a factor of the order of D. The diagram sequence represented in Fig. [TJis equivalent 



G 



to a closed system of nonlinear equations for the averaged magnetization (m;) and the 
correlation function Sfj a = (m a im a j) of the spin components with a = 3, .. . ,D. The 
diagrammatic equation in Fig. [l]a is the generalization of the Curie- Weiss equation 
(2.3) for the magnetization (the angle brackets are dropped): 



m, = dA(^, l m )/d£ t = A z (£i, l ai )e z + A y (^ i , l ai )e y 



(2.8) 



where the (averaged) molecular field £j is given by the expression (|2.2| ) without the 
last term and l a i is related to the dispersion of the molecular field fluctuations on the 
site i: 

lai = | = liVlP 2 E J H J i? S ?f- (2-9) 

jj' 

The spin cumulant averages A.., on a site i (see Fig. |l|a,c) renormalized by Gaussian 
fluctuations of the molecular field are given by 



D 



D 



a,/3=3 



\ .. A. f ^A,..„ Q l a + ^ ( 1-S a0 + -j£ a/3 ) A... mW U /9 + ..., (2.10) 



where taking into account only the first term corresponds to MFA. These series, 
describing the influence of pair-correlated fluctuations of different components of the 
molecular field, can be rewritten as 



A =11 E irrl^TSM- 



a— 3 n a — 



A...(£J = exp 



02 

E ^TjT? 

a=3 



A...(€i)- (2-11) 



Such exponential differential operators were considered by Horwitz and Callen |p8| for 
the Ising model. Generalization of their results for the multi-component case yields 
a closed formula 

A.. = ^/^ 2 re- 2 A...(0, (2.12) 
where Ci is the spread molecular field given by 



D 



(2.13) 



Q = 3 



and the integration in (2.12) is performed with respect to the (D — 2)-component 
vector variable r = {r a }. 

The Dyson equation for the spin-spin correlation function S?? entering fl2.9| ) is 



represented in Fig. [l]b and has the analytical form 



qua _ a r 



(2.14) 



Applying this equation two times, one can simplify the expression for the quantity l 
to 



In 



1 



2A C 



1 



(2.15) 



The system of SCGA equations for a domain wall in a biaxial ferromagnet stated above 
simplifies for a small anisotropy (1 — r\ <C 1) and in the limit _D — > oo, which will be 
pursued in the next section. For r\ a — the molecular field fluctuations measured by 
the quantity l a i vanish according to (2.14) and ( [2.15 ), and the magnetization equation 
(Ef) reduces to the MFA result (|2~3). 
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3. The spherical limit 



In the spherical limit D 
to 



the Langevin function (see (2.3), (|2.7|)) simplifies 



B(0 = 



1 



(3.1) 



where the scaled variable x (which should not be confused with the coordinate x 
appearing below) is given by x = Correspondingly, the spin cumulant averages 

(2.6) considered as functions of x can be estimated as 



1-fc 



(3.2) 



With the help of this estimate it can be shown that in the limit D — > oo SCGA 
becomes exact, since all other more complicated diagrams die out, at least as l/D 
|3jJ. Indeed, a unification of two oval blocks into a larger one, which leads to a more 
complicated diagram (e.g., A aa x App =>• A aa pp in Fig. |l]), leads to the appearance 
of an additional factor l/D, since 



(3.3) 



In pq | some of such higher-order diagrams were considered in the framework of the 
l/D expansion for low-dimensional classical ferro- and antiferromagnets. 

For the consideration of the limit D — > oo it is convenient to introduce the well 
behaved dimensionless temperature variable 8 = T/T^ FA , where J" c MFA = J /D 
and Jo is the zero component of the exchange interaction, as well as the further 
£>-independent variables: 



Gt = (D/9)A am , h = l a i/D, 
Expression ( ^.15 ) can now be rewritten as 



v = DS?f. 



20G l \8G 



(3.4) 



(3.5) 



and the expression for the square of the spread value of the argument x in (2.12) reads 

2 / \ 2 



2C 
D 



j 



Aii' 771 j 



2 >] 



16/ 



7 D 



Q = 3 



where Ay = Jij / Jo- It can be seen that the contributions of the fluctuations of the 
a-components of the molecular field to ( |3.6| ), each of them being small as l/D, are 
essential due to their large number. Now for D 1 the Gaussian integrals ( [2.12 ) can 
be easily calculated by applying the identity 



1-1/2 



f(ax 2 ) Si /(a/2), 



a < 1 



(3.7) 



successively D — 2 times. Thus, the integration leads simply to the replacement 



1/2 in (p.6[ ). Now for the quantity Gi of (3.4]with the use of the second of the 



formulae (2.7) and of the asymptotic expression (3.1) one gets 

1 



G* 



i + vTt: 



(3.8) 
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The magnetization equation (2.8) can be simplified by using the first of the formulae 
(|2~7|), as well as ([H]), to 



m z 



(3.9) 



Finally, determining x~ from (3.8) as a function of Gi, equating it to (|3.q) with 



1/2 and using (p.q) and (3.9), one arrives at the equation 



1. 



(3.10) 



which is nothing but the kinematic identity m? = 1 in the limit D — > oo. The 
normalized correlation function sw determined by (|3.4D satisfies the linear equation 



following from (2.14) 



(3.11) 



with the variable coefficient Gi. 



Equations (3.9)— ( 3.11 ) constitute the closed system of equations, which can be 
applied for the calculation of the domain wall magnetization profile in the spherical 
limit. 

In the homogeneous case (o r in o ne of the domains) m y — and m z and G are 
constants. In this case equation ( 3.11 ) can be easily solved with the help of the Fourier 
transformation, which results in 



Sii = i'O 



e?q 



s q = 9GP{r, a G), P(X) = v 



c/q 



1 



(2tt) 3 1-XXq 



(3.12) 



where vo is the unit cell volume and Aq = Jq/Jo- In the long- wavelength limit 
Aq = 1 — aq 2 , where a ~ a 2 , and a is the lattice spacing. The lattice integral P (X) 
has the following properties: 



n^J-i VF-c (l-X) 1 /2 ) l-X«l, 



(3.13) 



where z coincides with the number of the nearest neighbors for the nn interactions 
and W (the Watson integral) and cq = fo/(47ra 3 / 2 ) are lattice-dependent constants. 
For a simple cubic (sc) lattice vq = a 3 , and a = a§/6, hence cq = (2/7r)(3/2) 3 / 2 . 
Since the sum in the right-hand part of the first of the equations ([T^) equals m z , this 
equation is sati sfied only if m z — (above T c ) or G = 1 (below T c ). In these cases 
from equation ( 3.10 ) one gets the temperature-dependent magnetization m = m z : 



(1 



a/2 



< 



l/P{Va). 



(3.14) 



It can be seen that in the fully isotropic case {rj = r\ a = 1) the value of the phase 
transition temperature in the bulk 8 C = l/P(r) a ) reduces to the well-known result 

e c = i/w§. 

The width of a Bloch wall 5b in a uniaxial ferromagnet is determined by the 
balance between the anisotropy and inhomogeneous exchange energies. For small 
anisotropy (1 — 77 <C 1) the condition Sb a is fulfilled. In this case the change of 
the domain wall magnetization at a distance of ag is small, and for systems with nn 
interactions one can rewrite the lattice sum in (|3 .91) about some point r = as 



Aymj = m(r) + aAm(r), 



(3.15) 
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where the second term with the Laplace-operator A is small in comparison to the 
first one. The quantity Gi in ( |3.9| ) determined by (3.4) is in the small-anisotropy case 
also a continuous function of the coordinate with the scale S. Moreover, as can be 



seen from ( |3.9| ) and ( |3.15 ), the deviation of G from its bulk- value 1 is small, i.e., 

G(r) = 1 + Gi(r), d < 1. (3.16) 



Now one can rewrite the equations (3J3) in terms of the normalized m agnet ization 
n = m/m e , where m e is the equilibrium bulk magnetization given by (3.14). The 
result in the chosen geometry is of the form 



an"(i) = — Gi{x)n z {x), 

an"(x) = -Gx(x)n y (x) + (1 - rj)n y (x). 



The kinematic equation ( 3.10| ) can be represented in terms of n as 
e(l - n 2 ) = e cSii /9 - 1, e = 9 c /8 - 1. 



(3.17) 



(3.18) 



Unlik e the magnetization equation (3.9), the equation for the correlatio n fun ction 
( 3.11 ) cannot, in general, be written in a continuous form of the type ( 3.17 ). In 
the general case we are going to consider, where rj a is not necessarily close to 1, the 
correlation length of the a spin components below 9 C , 



£,ca = \/ar) a /(l - rj a ) 



(3.19) 



(not to be confused with components of the normalized molecular field £), which 
can be determined from sq in (3.12), can be comparable with the lattice spacing 
ap. Moreover, even in the case 1 — r] a -C 1 the continuous approximation for sw 



does not yield the correct bulk result (3.12) which is formed by the integration over 
the whole Brillouin zone and not only over the logwavelength region 9 « 1. But it 
can, nevertheless, be applied for the calculation of the wall properties, as we shall see 
below. 



4. Domain wall magnetization profile 



Before proceeding to the solution of equations ( |3.17| ), ( |3.18|) and ( 3.11 ) in the general 
situation, let us consider at first the case r\ a = 0, where the spin fluctuations play no 
role and the situation is described exac tly by MFA. Here the solution of ( |3.11 ) yields 
6Gi, and with the help of ( 0.18 ) one gets G\(x) = e[l — n 2 (x)]. Adopting it 



in (3.17) and solving the resulting equations, one arrives at the magnetization profile 
n z — ± tanh (x/S), n y = ±p/cosh(x/i5), (4-1) 

where 



P : 



VI - t, t = 2(1 -rf)/e< 1 
0, r>l, 



8 B = \ A*/(1 - rjj, 
5 L = ^Te = 5 B T 1 / 2 , 



T < 1 
T > 1. 

(4.2) 

It can be seen that the crossover from the Bloch (p = 1) to the linear (p — 0) wall 
proceeds with increasing temperature through the elliptic one having n 2 + n 2 / p 2 = 1 
with < p < 1, and the transverse component of the magnetization in the wall, 
n y , disappears through a second-order phase transition. The temperature of the DW 
phase transition can be written in three forms: 



t b =1, e B = 2(1 - 7]) «: 1, 



= l/(l + e B ) = l-e B , 



(4.3) 
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the quantity r playing the same role for a DW as the "absolute" temperature 9 for 
the bulk (cf. (3.14)). The temperature-dependent factor p in ( flip can be inter- 
preted as the DW order parameter m b |3l], Q ■ Whereas the Bloch-wall width 8b is 
temperature-independent, the width of the linear wall 8l is determined by the bal- 
ance of the homogeneous and inhomogeneous exchange energies and is diverging at 
9 C . Considering the first of equations (4.1) for x 3> 6l, one can identify 



1 - tanh(x/8 L ) = exp(-x/£ cz ) 



2Cc 



(4.4) 



where £ cz is the correlation length of the z spin components. One should also note the 
analogy between the Bloch-wall width 8b (4-2) and the transverse correlation length 
£ CQ (3.19), which coinside for a purely uniaxial (j] a — rj) model with small anisotropy. 
The function G\ entering equations (3.17) can be written as 



2a 



1 



cosh 2 (x/ 'S) S 2 cosh 2 (x/S)' 



(4.5) 



Since < 



in the small-anisotropy case G\ ~ (ao /8) 2 <C 1 in the whole temperature 



interval. 

Now we proceed to the solution of the magnetic interface problem described by 
equations ( [3.17 ), Q3.18 ), and (3.11) in the general case rj a ^ 0. The solution of the 
Dyson equation for the correlation function (3.11) depends on the relation between 
the correlation length £ CQ of (3.19) and the other length scales, ao and 8. If £ CQ 3> a , 
which is satisfied for 1 — r\ a <C 1, the continuous approximation to equation ( |3.11 ) 
can be applied. If £ CQ <C 8 (i.e., 1 — rj a ^> 1 — rj, see fl4.2|)), the correlation function 



Sij can be easily calculated locally with respect to the slowly changing magnetization 
profile (or, more exactly, the profile of G) in the wall. For uniaxial ferromagnets with 
a small anisotropy (1 — r) <C 1) considered throughout this paper these limiting cases 
overlap in the region \ — rj^il — rj a <^il. Let us con sider at first the case £ CQ <C 8. 
Here one can make a replacement Gi => Gy in (3.11), after which this equation can 
be solved like in the homogeneous case. With the use of |[l|) and d3T6|) one gets 



£* OGiP^Gi) = 9P( Va )[l + J(i, )Gi(x)], 



where 



r/gP' {Va) 
Pfrcc) 



1 + 2?? 2 /z, 
CQ 1 

2P(Va) VT^fo' 



i - na < i 



(4.6) 



(4.7) 



Gi through the magnetization profile n(x): 

G 1 (x)=e[l-n 2 (x)}r 1 (r ]a ). 



and Co is determined by (3.13). Now with the use of (4.6) and (3.18) one can express 

(4.8) 



This expression differs only by a constant from that used in the beginning of this 
section in the MFA limit r\ a = 0. Solving now the magnetization equations ( |3.17| ) 
like above, one gets the same DW magnetization profile ( |4.1| ), where the parameters 
p and 8 are given by (4.2) with the renormalized DW temperature: r => f = rJ(r? a ). 
The critical values of the three temperature variables (cf. Q4.3| )) read now 

T B =I- 1 (r ]a )<l, e B = 2(1 - n)Hn*), 8b =0 c /(l + e B ), (4.9) 



where 9 C is given by ( |3 .14 ). One can see from (4.9) and (4.7) that for 1 — r\ a <C 
1 the effective transition temperature tb becomes small. On the other hand, due 
to the validity condition 1 — r\ a ^> 1 — r\ the absolute temperature 9b of the DW 
phase transition remains in the limiting case under consideration close to 9 C (e^ <C 
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1). It can be seen that in this case the domain wall does not demonstrate any 2- 
dimensional behaviour, and its phase transition at 9 b < 9 C can still be described 
qualitatively by the effective diminishing of the ordering interaction for the wall spins 
forced perpendicularly to the easy axis z, as was said in the Introduction. The effect 
in the case 1 — r\ a <C 1 is much larger than according to the MFA estimates because 
of the divergence of the function I(rj a ) of (4.7). Accordingly, the width of the linear 
wall can considerably exceed its mean-field value (4.2): 




P(Va) (1 - Va) 1 ^ 



(4.10) 



in the case 1 — rj a -C 1. Since 8l is related to the longitudinal correlation length 



(see (1.4)), this result shows a non-trivial influence of fluctuations of transverse spin 
components on the longitudinal spin correlations in the anisotropic spherical model. 
One can also check that the function G\{x) of (4.S) is still given by the expression 
(4.5) with the changed value of the DW width 8. 

In the other limiting case, £ CQ 3> eto, a continuous approximation of the type ( [3.15 ) 
can be applied to the Dyson equation ( [3.11 ) . Making the Fourier transformation with 
respect to the coordinates y and z and using the conditions 1 — r\ a <C 1 and ( [3 . 1 ) 
one arrives at the differential equation for the correlation function 



as" (a;) - [1 - r) a + aq 2 - G!(x)]s(x) = -9S(x - x'), 



(4.11) 



where q 2 



q 2 and the "mute" argument x' of s was dropped. This equation 



should be solved to yield s with x — x 1 as a function or functional of Gi 1 and Sa 
(cf. (|4.6|) ) should be obtai ned b y the integration of s over q y and q z . Then, as a bove, 
G\ should be found from ( 3.1 8| ) and used in the magn etization equations ( 3.17 ). All 



this seems to be too complicated since equation (4.11) cannot be solved analytically 
for the arbitrary function G\(x). But the expected result that the DW transition 
temperature 9 b turns to zero in the purely uniaxial case r\ a = 77 signals that there 
should be an exact solution to the problem. We can try to find it assuming that 
G\(x) has the same functional form as above, (4.5), with some renormalized value of 
the DW width 8 as a parameter. Then using a new variable u = tanh(x/(5), equation 
(4.11) can be rewritten as 



Ul-u 2 )^- 
du du 



2- 



/< 



1 



s(u) — 5(u— u'), 

a 



I* 



and solved in terms of the adjoined Legendre functions 



which leads to 
f 

s(x,x',q) = 



2au 



exp 



(-fk-*1 



u ± u f 1 + u 



T{2±n) \l-u 



— {l-^+aq 2 ) (4.12) 
a 



(4.13) 



1 - tanh § tanh % / \x — x'\ 

1 H I 1 + H tanh 1 



u 2 - 1 V ' ' 6 

(4.14) 

This (not translationally invariant) expression can be reduced in the ' with 

the help of (|4.5|) to the form 



s(x,x,q) 



2a 1 / 2 v / l ~rj a + aq 2 



Gi(x) 



2[1 - Va -(1- rj)/5 2 + aq 2 } 



(4.15) 



where 5 = 5 /8b (sec (4.2)). Since the structure of this expression is analogous to that 
of (4.6), it is now clear that the choice of G\ in the form (4.5) was right. Integrating 
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(4.15) over the 2-dimensional wavevector q to get su and proceeding as above, one gets 
functionally the same results (4.1) and (4.2) with a new renormalized DW temperature 



= Tl(r] a ,r),6), 



where 



CO 



(4.16) 



(4.17) 



The latter simplifies in the l imit 1 — rj a 3> 1 — ij to the second limiting expression in 
(4.7). It can be s een from ( 4.15 ) that in the case under consideration, 1 — r\ a <C 1, 
the integral ( 4.17 ) is determined by the long- wavelength region, which justifies using 
the continuous approximation for the transverse correlati on fu nction s. In the region 
of elliptic walls (f < 1) one has 5 = 5 b and hence in ( [4.17 ) one has 5 = 1. The 
critical values of both the DW temperature r and the absolute temperature 9 ( |4.9| ) 
determined now by I(rj a , f], 1) go to zero in the uniaxial limit: 



tb oc V b 



oc 1 / In • 



1 



Vet -> T)> 



(4.18) 



which corresponds to the 2-dimensional nature of domain walls. One can see that the 
2-dimensional effects lead to a further decrease of the DW phase transition tempera- 
ture: 



t b =I 1 {Va,V, 1) < t b } = 1 1 ( 1 la), where I(r] a ) is given by (O) 



In the linear- wall region (f > 1) the normalized wall width 5 is given by the 
solution of the transcendental equation 5 = Tl(rj a ,r),S) following from (4.2) and 
(|41§). The latter can be rewritten in the natural units in the form 



5l = 



co 



In 



5l+S,c 



2e P( Va ) 5 L - e 



(4.19) 



where £ cct is the transverse corr elation length given by ( 3. IS ) . Far from tb , where 
Sl 3> £ca, the solution of ( 4. IS ) leads to the formula ( 4.10 ). This asymptotic de- 
pendence can also be represented in the form 5 = ^jTl(rj a ) = \J t/t b °\ which is the 
analogue of Curie- Weiss asymptote for the bulk susceptibility x{T) of a ferromagnet 
far above T c . In the purely isotropic model, rj a — rj 1 the linear DW structure is real- 
ized in the whole temperature range. The DW width 5l determined by the solution 



of (4.19) shows a crossover to the Bloch-wall width 5b at low temperatures: 



5l=5 l 



1 + 2 exp 



2e P(r/) 



a/1 - rj c 



e = Oc/0 - 1 > a/1-77. 



(4.20) 



It is worth to note that in this limiting case the longitudinal correlation length, £ cz = 
6l/2, is determined by the transverse one, £ CQ = 6b- 

The temperature dependences of the DW order parameter = p and the 
normalized DW width 5 are represented in Fig. for different values of ?7 Q and 
1 — i] = 10~ 3 . One can see that for 1 — ?7 a = 10~ 2 the fluctuational decrease of 
the DW transition temperature is essential, although 2-dimensional effects are still 
negligible in this case. In contrast, for 1 — r\ a = 1.2 • 10 -3 they come into play, the 
corresponding additional diminishing of the transition temperature becomes essential, 
and the temperature dependence of 5 in the logarithmic scale is no longer a straight 
line. 

5. Discussion 

In the main part of this paper the structure of domain walls in the biaxial ferromag- 



netic model described by the classical- vector Hamiltonian (2.1) was calculated exactly 
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Fig. 2: Temperature dependences of the domain wall order parameter raj = p and 
the DW width 5 for different values of the in-plane anisotropy parameter r\ a . 



in the whole temperature range T < T c in the spherical limit D — > oo. The main 
qualitative result is that in the purely un iaxia l model (rj a = rf) the Bloch walls do not 
exist at any non-zero temperature (see ( 4.18 )) being disordered by thermal fluctua- 
tions to the linear (Ising-like) ones. This result complements the well-known fact that 
the Bloch walls in the purely uniaxial model (considered within the phenomcnological 
micromagnetic approach which is equivalent to MFA) cannot move since their maxi- 
mal velocity (the Walker velocity) is equal to zero. In the opposite limit, r] a = 0, the 
model with D — > oo total spin components and a finite number (here 2) of interacting 
ones is realized. In this case spin fluctuations die out and MFA becomes exact. The 
temperature of the phase transition from Bloch to linear walls Tg changes as function 
of rj a from its MFA value T B = (1 - e B )T c , e B — 2(1 - n) < 1 at r] a = to at 
f] a = r) (see Q and Q). 

In the actual case of a small anisotropy, 1 — rj <C 1, the behavior of a domain 
wall is more complicated than that of a purely 2-dimensional object, since the DW 
width S is much larger than the lattice spacing ao- In the spherical limit this leads 
to the existence of two different mechanisms of the DW ordering, depending on the 
value of n a . For the strong anisotropy in the basis plane x, y, i.e., 1 — r\ a 1 — r), 
the correlation length £ CQ of (3.19) for all temperatures is much shorter than the wall 



width 6 in ( |4.2| ) , and the DW phase transition at Tb can be inte rp rete d as the locally 
shifted bulk one (the "perturbed 3-dimensional" case, see 0)-©). In contrast, 
in the case of the two anisotropies comparabl e wit h each othe r, 1 — r\ a ~ 1 — rj, the 
true 2-dimensional situation is realized (see ( 4.17 ) and ( 4.18 )). Such a separation 
does not, however, persist for models with finite values of D (e.g., for the Heisenberg 
model, D = 3), where the DW phase transition should always be of a 2-dimensional 
character. For such models the fluctuations of the ordering spin components, m z and 
m y , also play a role, and in the temperature interval about Tb, where the diverging 
correlation length £ cy exceeds the DW width 5, a 2-dimensional behavior is realized. 
Since 8 3> ao, this temperature interval should be much narrower than for a purely 
2-dimensional system. The asymptotic critical behavior of the DW order parameter 
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tub — p in (4.1) is described by the critical index (3b = 1/8 of the 2-dimensional Ising 
model, as was confirmed experimentally in |3l], |3^] . 

As a subject of future investigations, the temperature dependence of the DW 
magnetization profile in the self-consistent Gaussian approximation, without going 
to the limit D — > oo, can be considered. Although it can be connected with more 
complicated numerical calculations, one can expect to obtain, with the help of SCGA, 
essentially more accurate results for domain walls in comparison with the spherical 



approximation, as was demonstrated for the bulk properties [|40, 18 . The other prob- 
lem is to formulate dynamic equations for fluctuating domain walls and to calculate 
their mobility. Such equations can, in principle, be obtained with the help of some 
dynamical generalization of the classical spin diagram technique ]j~7|| . Unfortunately, 
in the dynamical case one cannot make use of going to the limit D — > oo with all 
related simplifications, and only SCGA for the Heisenberg model can be used as the 
underlying static approach. 

A promising field for the application of the methods of this paper is the surface 
effects in finite and semi-infinite magnetic systems, which are very sensitive to the 
anisotropy. This problem was addressed recently in [fi2[ . 
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